Coulomb breakup of ^Li into a + d in the field of a ^°^Pb ion 
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The triple differential cross section of the ^''*Pb(®Li, Q:d)^''*Pb quasielastic breakup is calculated 
at a collision energy of 156 MeV and a scattering angle range of 2° -6°. We fit the parameters of 
the Woods-Saxon potential using the experimental a-d phase shifts for different states to describe 
the relative motion of a particle and deuteron. To check the validity of the two particle approach 
for the a-d system, we apply a potential model to describe the ^H[(a,7)^Li radiative capture. We 
calculate the Coulomb breakup using the semiclassical method while an estimation of the nuclear 
breakup is made on the basis of the diffraction theory. A comparison of our calculation with the 
experimental data of Kiener et al. [ Phys. Rev. C 44, 2195 (1991)] gives evidence for the dominance 
of the Coulomb dissociation mechanism and the contribution of nuclear distortion, but is essentially 
smaller than the value reported by Hammache et al. [ Phys. Rev. C 82, 065803 (2010)]. The 
, results of our calculation for the triple differential cross sections (contributed by the Coulomb and 

■ nuclear mechanisms) of the ®Li breakup hint toward a forward-backward asymmetry in the relative 

direction of the a-particle and deuteron emission, especially at smaller scattering angles, in the ^Li 
center-of-mass (cm.) system. 
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I. INTRODUCTION 



O , The study of processes, relevant for nuclear astrophysics, by indirect methods gives the possibility to extract 
^ ' the astrophysical S factor at extremely low energies when extraction by direct methods is not possible due to the 
, ^ , 1 Coulomb barrier suppression. Among these indirect methods we cite the elastic Coulomb dissociation method sug- 
gested by Baur, Bertulani and Rebel iT-'s'], the asymptotic normalization coefficient method (ANC) suggested by 
fT^ \ Mukhamedzhanov and Timofeyuk [4-ii], and the Trojan horse method (THM) suggested by Baur and modified by 
K*" ■ Spitaleri 0, The study of nuclear reactions at high energy is, in general, very complicated owing to the strong 
nuclear interaction between the colliding nuclei. However, in the peripheral collisions of a light nucleus with heavy 
target the reaction mechanism becomes simple owing to the negligible contribution of nuclear distortion and excitation 
becomes purely Coulombic. Electromagnetic excitation is a very powerful tool for the extraction of the information 
concerning radiative capture at extremely low energies when the direct measurements of the radiative capture of 
1^ ■ nuclei for astrophysical purposes is impossible. The Coulomb dissociation ex per iments are being performed at dif- 
2 fcrent centers around the world (e.g. GSI, Germany [9j; MSU/NSCL, USA [TO]; RIKEN, Japan The value 

I ■ of the astrophysical S factor extracted from the Coulomb breakup experiment is affected by various uncertainties, 
• • ' namely: (i) the method of extrapolation of the data to zero energy; (ii) contributions from various electromagnetic 
. ^ multipoles (E'l, Ml); (iii) assumptions about the nuclear interactions; and, (iv) various higher-order effects (see 
' Refs. |13 - [la ] and references therein). The astrophysical S factor at zero energy extracted from the direct radiative 
5^ , capture and the Coulomb breakup reaction should be the same. However, the extrapolation to zero energy from the 
_<^: d{a, -ffU and 208pb(6Li, a d)^^^Ph reactions gives different values for the astrophysical factor. This is owing to the 
fact that at ad relative energy below 100 keV the El dipole cross section becomes larger than the E2 quadrupole 
cross section for the direct capture process, while for the Coulomb breakup the El cross section is always smaller than 
the E2 one at any given energy. Kiener et al. [l^ investigated the Coulomb breakup of ^Li in the field of ^°^Pb ion 
and extracted large value of the astrophysical factor at zero energy from extrapolation in the range Ead < 400 keV. 
Shyam et al. [l7j underlined that the result of the authors of Ref. il6| was free from the nuclear background which 
maybe important. Kiener et al. did not take the contribution of the dipole transition to the cross section owing to 
isospin selection rule. However, this rule is violated and El transition can still occur. Recently a new measurement of 
the Coulomb breakup at extremely high energy of ^Li (150 A MeV) was performed [Ts] and the authors have claimed 
disclosure of evidence for the large contribution of the Coulomb-nuclear interference. However, the authors did not 
present the cross sections (the histograms for counts in Figs. 6 and 7 of Ref. [IQ] were shown). Figure 10 of Ref. 
[18| showed the ratio of nuclear to Coulomb differential cross sections for ^Li which is very large compared to the 
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qualitative estimation given earlier in Ref. [19| . 

The dissociation by the nuclear field of a target can be excluded by observation of the fragments of the reaction at 
forward scattering angles. This extremely small angle corresponds to the impact parameter essentially larger than the 
sum of the nuclear radii of the projectile and the target. For instance, in the experiments fl^ the angle was varied in 
the range 2°-6° corresponding to the impact parameter in the range 20-65 fm. However, the dipole transition in the 
^Li Coulomb breakup is suppressed as in the d{a, 7)^Li radiative capture, therefore, we can expect that the nuclear 
dissociation gives a relatively large contribution to the ^Li Coulomb breakup cross section. 

The problem of the higher-order effects was discussed in many papers (see, e.g. Refs. [1, [l^], and references 
therein) . In the semiclassical theory these effects are inversely proportional to the impact parameter and the velocity 
of collision. Therefore we can constrain ourselves by the first-order amplitude for the excitation of a fast nucleus at 
sufficiently small scattering angles. Among higher-order effects, we also note the three-body Coulomb effects in the 
final state, which were discussed in Refs. |2ll |23]. 

In this paper we consider the dissociation ^°®Pb(^Li, a d)^°*Pb using the time-dependent perturbation theory for 
the Coulomb breakup, whereas the nuclear breakup is viewed as a diffractive dissociation. Our treatment is used 
for the complete analysis of the experiments, which were performed with a 156 MeV ^Li beam at the Karlsruhe 
Isochronous Cyclotron [16| . This reaction is relevant to the c?(a, 7)®Li radiative capture, which is one of the important 
nucleosynthesis reactions. When the collision energy is sufficiently large, we can consider the motion of the center of 
mass of a projectile along a classical trajectory. At small scattering angles of the projectile the value of the impact 
parameter would be so large that at such distances only the Coulomb interaction between the projectile and the target 
would be significant. The nuclear breakup does take place, but it mainly occurs near the target surface where the 
nuclear potential of the interaction falls down quickly. At small scattering angles and high energy collisions, the most 
appropriate method of estimation of the nuclear breakup is the diffraction theory (23| as also has been mentioned 
above. Here we mention the experiments concerning the ^Li breakup at energies 60 MeV [l^l and 31, 33, 35 and 39 



MeV |25|, which also confirmed the Coulomb dissociation mechanism of Li on ""Pb. However, from our opinion the 



authors of Ref. [24| made an erroneous conclusion that the second-order Coulomb excitation theory can improve the 
deviation in the angular distribution between the theory and the experiment. On the basis of the developed method 
we desire to make the analysis of the data presented in Ref. [l^ and to explore the possibility of the extraction of the 
astrophysical S factor from the ^''*Pb(^Li, a d)^°^Pb reaction by a comparison of our results with the experimental 
data of Ref. [l|. 

We use the system of units in which h — c — 1. 



II. GENERAL FORMALISM 



In the semiclassical theory the center-of-mass (cm.) motion of the projectile is considered classically, whereas 
the relative motion of the clusters in the projectile is treated completely in a quantum mechanical fashion Q. In 
this approach the purely Coulomb breakup cross section of peripheral ^°^Pb(^Li, Q;d)'^"^Pb reaction is presented as a 
product of the Rutherford scattering cross section of ®Li in the field of ^°^Pb ion and the probability of the ®Li -^a + d 
disintegration. We calculate the probability of disintegration in the ^Li frame (projectile frame) because the calculation 
does not depend on any reference frame. In this frame a heavy ^°^Pb moves along a straight line with a constant 
velocity v. If the energy of the ad relative motion is small we can restrict by using only the low partial waves at 
I = 0,1,2 for the description of the ad motion. The E2 multipole gives the main contribution to the transition 
amplitude of the electromagnetic dissociation of the ^Li, however, the El transition should also be included to the 
amplitude owing to a violation of the isospin forbidden rule. 

The time-dependent perturbation for the Coulomb breakup A + a^A + c + b is 

where Zag is the charge of a heavy ion A (^°^Pb target), the projectile a (^Li nucleus) is dissociated into b (deuteron) 
and c {a particle), R(t) = h + vt gives the position of the target in the projectile frame, and p(x) is the charge 
density operator, b is the impact parameter. For the peripheral reaction we can take the charge density operator in 
the two-body approach as 

p(x) = Zbe(5(x - rfc) + Z^eSi^ - rj, (2) 

where defines the position of particle i in the projectile frame, and ZiC is its charge. 

In perturbation theory the amplitude of the transition from the initial state \i > (wave function of a = b + c system 
in the ground state) to the final state \f > (wave function of the b + c system in the continuum state) is given as a 
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«/. = Sf, + ay> + ay> + • • • , (3) 

of different order contributions. We use the first-order constrain for the amphtude of the transition. Expanding H{t) 
in muhipoles we get 



ji) _ iirZAe ^ (-1)^ 
where M{X, fi) is the electric multipole operator 



A^(A,/i) = / d^xp(x)x^yA^(x). (5) 
The semiclassical orbital S\^{uj) integral is given by 
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where x and tt(t) are the unit vectors along the position vectors x and R(0. respectively. The integrals Sx|^^{i^J) were 
calculated analytically and the results of calculations were presented in Ref. [2a] . Determining the relative coordinate 
r = Tb — Tc, inserting Eq.® into Eg.© and performing the integration over variable x we obtain 



M{X,^,)^^^t ^ + (-ir^ r^Yx.ir), (7) 

where fibc is the reduced mass of particles b and c. We see that the dipole transition operator /i) is not equal to 

zero because the value of the charge-to-mass ratio is slightly different for the a particle and deuteron. Therefore, the 
El transition amplitude is not zero. 

The triple cross section of the Coulomb breakup from the ground state with angular momentum and parity = 1+ 
of ^Li to the final state with relative momentum k of the a particle and deuteron having reduced mass ^ad can be 
expressed in terms of the excitation amplitude a fi : 

(Pac dan 1 i „ 12 ^^adk 



driaddnudEc^d d^Li 2 J, + 1 ^ ' ' ' (27r) 

The elastic Coulomb cross section da-ji/dflu is calculated classically for the scattering of the cm. of the projectile 
6Li. 

Next we consider the nuclear breakup. As mentioned earlier, the nuclear breakup occurs mostly near the surface 
of a heavy target nucleus ^'^^Pb. If the scattering angle is small we may apply the diffraction theory assuming that 
the target nucleus is a completely absorptive "black" sphere of radius Rbi- The amplitude for the elastic breakup 
according to the diffraction theory is given by 

^(q, k) = ^ 1 d^e^^'' J d\r^{rMh, r)^o(r), (9) 

where q = — k/, is the initial momentum of ^Li and k/ is the final momentum of the center of mass of the ad 
system after dissociation; | i) — '!/'o(r) and | /) = '0k(r) are the wave functions of arf in the bound and continuum 
states, respectively; b is the impact parameter of the ad system; uj{h,r) is the total profile function for the ad 
system. We may take the vector q to be orthogonal to the momentum vector k; due to high energy collision and 
small scattering angle (g « fc^ • 9, where 9 is the scattering angle). The total profile function 

CJ Wq + W<J - WaCJrf (10) 

is composed of the profile functions of the fragments. The third term in Eq. (fTO| describes the double scattering 
and its contribution to the cross section is much smaller than the contributions of the first two terms. Neglecting the 
double-scattering term we obtain the cross section of the nuclear breakup as 

fcji?bi s(q, k)| 7^-^, (11) 



dfladd^LidEad q (27r 
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where Ji (x) is the Bessel function and 



s(q,k)= j d're^^''r^{r)Mr)- (12) 

The total cross section is considered as the sum 

d^crt _ (Pac ^ d^CFN ^^^^ 

dQaddflhidEad dfladdiludEad dfladdfludEad 

In Eq. ([T^ the Coulomb-nuclear interference term is neglected. We note that the Coulomb and nuclear breakup are 
calculated using different approaches therefore we can not calculate the interference term. Even if both the Coulomb 
and nuclear breakup amplitudes were calculated using the same approach, one expects the interference term to be 
small (see for example the calculation in the distorted wave Born approximation (DWBA) method used by Bertulani 
and Hussein, Figs. 4 and 5 in Ref. [13 )■ We further note that at the small scattering angles considered in this 
calculation, the elastic Coulomb scattering cross section calculated classically is equal to the cross section calculated 
by the diffraction theory. 



III. a-d POTENTIALS 



To describe the relative motion of the a particle and deuteron we use the Woods-Saxon potential with the orbital 
terms, 
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VNir) = -Vo- Vsiil ■ s)^- — ^ (14) 

mt r ar u _ _ w „ i 



1 -I- exp [(r - RN)/a] ' 



with the standard value of the diffuseness a — 0.65 fm. The parameters Vq and Vgi for different ad states is fitted 
from the shift phase analysis. We take the nuclear radius of the potential Rn as 

i?Ar = ro-Ai/3, (15) 

with the standard value tq = 1.25 fm, and the nuclear mass number A = 6 for ®Li . The Coulomb potential is taken 
as 

I r>Rc, 

where ZaS and Z^e are the charges of the a particle and deuteron, respectively; Rc = rc A^^^ {rc = 1-25 fm). 

The parameters of the depth of the potentials Vb and Vsi were fixed by fitting the experimental 5", P and D phase 
shifts of the elastic a-d scattering [28l - l32j and the binding energy of the ^Li ground state. We obtained the following 
values of the depths: Vb = 60.73 MeV for the state; Vq = 57.0 MeV, V^i = 4.0 MeV for the ^Pi, ^Pi and "^Pi 
states; Vb = 55.9 MeV, V^; = 4.0 MeV for the ^Di, "^Di states and Vb = 55.9 MeV, Vsi = 5.06 MeV for the ^Di 
state. To describe the 3^ resonance of ^Li correctly, the depth of the spin-orbital part of the potential for ^Di state 
is taken slightly differently. We note that our fitted parameters of the Woods-Saxon potential are slightly different 
from the parameters used in Ref. [l^. This difference is obvious as the phase shifts are determined with errors and 
the resulting fit can give rise to such differences in the parameters of the potential. Note the spin-orbital potential in 
Eq. ([T4| contains the dimensional parameter 1/mJ = 2.136 fm^ while the same potential in Ref. [l^ has = 4 fm^. 

It is important to remember that at low energies the ad radiative capture depends on the tail of the ^Li bound 
state wave function projected on the a-d channel [ssl |33 | . The amplitude of this tail is the asymptotic normalization 
coefficient (ANC). We note that the range for the ANC obtained by various techniques is wide (C = 1.51-3.25 fm^^^^) 
[35} . We calculate the amplitude of the peripheral radiative capture reaction ^II(a, 7)^Li at very low collision energies 
where the ANC for the virtual decay ^Li a -\- d from the ground state governs the overall normalization of the 
peripheral reaction cross section. The value of the ANC (Co) obtained here using the fitted potential for the state 
of ^Li is 2.7 fm~^^^. This value is larger than C = 2.3 fm~^^^ extracted from the elastic a-d ^Si experimental phase 
shift by the analytic extrapolation to the pole of the partial scattering amplitude corresponding to the ^Li ground 
state [3^1. The same value of the ANC was also obtained from the solution of the three body a-p-n equation j36| , 
which was used in Refs. [H, [s^l and confirmed recently by ah initio calculations [33|. To obtain the ANC value of 2.3 
fm ' we can find the phase-equivalent potential and the corresponding wave function by the method described in 
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Ref. Ref. [38j and discussed in Ref. [33|. The phase-equivalent potential does not change the scattering phase shift 
and the binding energy, but allows one to get the needed value of ANC. The phase-equivalent potential has the form 



where fcsir) means 



fes{r) = In 



1 + (A-1) 1- / u'^{r)dr 



The corresponding new wave function of the bound state is equal to 

u{r) 



Mr) = X'/' 



l + (A-l)/>2(^)d^' 



(17) 



(18) 



(19) 



where u(r) is the wave function obtained from the solution of the Schrodinger equation with the parameters of the 
Woods-Saxon potential determined from the phase-shift analysis. We note that, to decrease the value of the ANC 
we have to take A > 1. For our phase-equivalent potential and the wave function corresponding to C = 2.3 fm~^/^ 
the value of A is 1.38. [This follows from Eq. (|19p : C = A~^/^Co]. A detailed description of this method and its 
application in case of the a-d radiative capture can be found in Ref. [1^. Thus for calculation of the cross sections 
of the ^H(a,7)^Li radiative capture and the ^*'^Pb(®Li, a d)^°*Pb breakup we use the Ucs{r) bound wave function 
having C = 2.3 iiR-^^^. 



IV. RADIATIVE CAPTURE REACTION ^R{a, jfU 

Experimental measurements of the cross section of the ^H(a, 7)®Li reaction at extremely low energies are very 
difficult to carry out because the cross section is the order of a few nanobarns and decreases exponentially if energy 
goes to zero. The experimental results for the cross sections of the direct d + a capture were measured by Robertson 
et al. jsl] at cm. energy range 1-3.5 MeV, by Mohr et al. [i^l at the resonance point of 711 keV, and Cecil et al. 
(ilj at an ad cm. energy of 53 keV. Furthermore, we mention Refs. [i^ HI] where the analysis of the experimental 
results and theoretical calculations are given. Nollett et al. (i^l used a six-body approach for the calculation of the 
ad capture. From the analysis of the results at the E^ii energy close to zero we see an essential difference in the value 
of the astrophysical S factor depending on the applied value of the ANC. At extremely low energies the initial wave 
function of the ad system ceases to depend on the parameters of the nuclear potential, which is used to calculate this 
wave function, since the nuclear scattering phase shifts tend to zero. Therefore, we can replace the wave function of 
the initial state of the purely Coulomb wave if E^d < 100 keV. For any radiative capture at low energy, the main 
contribution comes from El, E2 and All transitions. However, for the ^H(a,7)^Li reaction the main contribution 
comes from the E2 quadrupole transition at the energy larger than 200 keV. The El transition begins to dominate if 
the Ead energy becomes less than 100 keV. The Ml capture remains negligible for all astrophysical interesting range 
of energies. 

For the calculation of the d{a, 7)^Li cross section we use the Woods-Saxon potential with the parameters described 
in the previous section. The initial wave function of the ad system includes the P and D waves which are solutions of 
the Schrodinger equation. It is clear that the ad direct capture is a peripheral reaction at the low energy (E^d < 300 
keV). Accordingly, the cross section is not sensitive to the choice of the parameters of potential describing the 
continuum states. It is rather strongly dependent on the ANC of the ground state. Such a property was used for 
the ANC calculation [sl, [33| , where the asymptotic wave function of the ^Li ground state wave function in two-body 
approximation was applied to find the astrophysical S factor. The main contribution to the matrix elements of the 
direct reaction at E^d < 300 keV comes from the external part of the used wave functions, while the internal part 
gives a very small contribution. Figure [T](a) shows the radial part of the integrand for transition from to '^Si at 
the energy Ead —0-1 MeV. We can see that the replacement of the P continuum wave function by the corresponding 
regular Coulomb wave function leads to almost the same value of the matrix element for the '^Pi -^^ Si transition. 
We further note that the same conclusion can be made for the transition from D states at the low energy. 

It seems that in the resonance region the complete microscopic model should be used (see Ref. [i^ and references 
therein) because the contribution to the transition amplitude from the internal part of the radial wave functions is 
also expected. Nevertheless, the two-particle approach can be used to describe the ^H(a,7)^Li resonance reaction 
and to get a quantitative result for the cross section at the resonance energy region. Such instances maybe explained 
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FIG. 1: (Color online) (a) The integrand of the matrix element for transition from "Pi to ^5*1 at energy Ead = 0.1 MeV 
calculated with the continuum wave function which is the solution of the Schrodinger equation (solid curve) and with the 
regular Coulomb wave function (dashed curve), (b) Same as in (a) for the transition from ^_Di to ^Si at the resonance energy 
E^a = 0.711 MeV. 

by the suppression of the El transition. The resonant amplitude is the result of the transition from the state '^Di to 
the state '^S'l and the transition remains peripheral due to the large centrifugal barrier in the ^Di state. The result 
of our calculations shows that the two-body wave function for the ground state of ^Li gives an acceptable result for 
the resonance cross section and the internal part of the wave function gives very small contribution to the transition 
amplitude at the resonance energy [see Fig. [D (6)]. The result of the calculations of the astrophysical S factor as a 
function of energy is shown in Fig. [5]over a wider energy range than shown in Fig. 3 of Ref. |33j] . It is clear from Fig. [5] 
that the two-body approach for the description of the ground and continuum states of the ad system gives perfectly 
good results at the low energy, including the 3"^ resonance energy region. The value of our calculated astrophysical 
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FIG. 2: (Color online) The calculated astrophysical S factor compared with the results of other theoretical calculations and 
experimental data. The solid line is our result for the total S factor, the dashed and dot-dashed lines are the contribution 
of E2 and El multipolarities, respectively. The dotted line shows the result of the calculation of the authors Ref. and 
dot-dot-dashed line corresponds to Ref. [1^. The experimental data (A, •, x) were taken from their graphical presentation in 
Ref. [H. 

S factor coincides with the one presented in Ref. [s^] for the low energies. At an energy larger than 1 MeV our 
result, definitely, agrees well with the experimental data [s^, while the results of the authors of Refs. [T8l.l43l| clearly 
overestimated the data. We also see a big disagreement of our result with the astrophysical factor extracted from the 
Coulomb breakup experiment [l^ below 500 keV. We note that the ANC determines the value of the astrophysical 
S factor at the energy range Ead < 300 keV where the reaction has a clear peripheral mechanism. The ratio of the 
astrophysical S factors (or cross sections) calculated near zero E^d energy using the ground wave functions with a 
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different value of the ANC is equal to the square of the ratio of the corresponding ANCs. As was also mentioned 
above, the ground state wave function of ^Li used in our calculation for the radiative capture has an ANC equal to 
2.3 fni^i/2 ^j^iig ^j^g ANC of the bound state wave function used in Refs. [Uliil is 2.7 fm-^/^ ^^^^ 3 2 
respectively The secondary peak of the calculated cross section near Ead = 3 MeV appears due to the wide 

resonance in the scattering wave at Ead — 2.838 MeV. The applied potential describes this resonance. To explain 
the disagreement of our results with the experimental data for Ead > 3 MeV, we refer to the work of NoUett et al. 
[4^, where the nature of this disagreement is discussed in detail. Additionally, we note that for the energy Ead > 3 
MeV the reaction is no longer peripheral. We have not calculated the cross section above 6 MeV because the phase 
shifts have a large imaginary part due to the open a + p + n, ^He + p and ^Li + n channels from the energy E > 5 
MeV and the restriction by the single ad channel becomes incorrect. From Fig. [2] we see that the El cross section 
dominates the E2 cross section at energies below 100 keV and Sei{0) = 1-01 MeVnb while Se2{0) = 0.21 MeVnb. 
Hence our total calculated value of the astrophysical factor at zero energy is 5(0) = 1.22 MeVnb. 

V. RESULTS FOR THE COULOMB BREAKUP OF ''Li 

Using the El and E2 multipole matrix elements calculated to determine the ad capture cross section, we compute 
the ^°®Pb(^Li, a d)^°^Pb breakup reaction. The calculated triple differential cross section [Eq.®] includes El, E2 
and E1E2 terms. The analyzed experimental data of the Coulomb breakup were taken from Rcf. [16], where data are 
presented for the scattering angles Quh — 2, 3, 4, and 6° in the laboratory frame in Tables III theough VI; and Figs. 7 
and 10, respectively. To convert the laboratory cross sections into the ad cm. cross sections we used transformation 
method described in Ref. [31 ■ Previously, only the data at the scattering angle 8iab = 3° were analyzed and the 
astrophysical factor at zero energy was extracted from that analysis [l6l |. The experimental data at other scattering 
angles have not been analyzed so far to the best of our knowledge. We chose the impact parameter corresponding to 
the selected scattering angles when we calculate the Coulomb breakup contribution. In spite of the fact that the El 
cross section becomes larger than the E2 cross section for the direct radiative capture at energies Ead < 100 keV, 
the El triple cross section is always less than the E2 one for the Coulomb breakup at any energy including the range 
Ead < 100 keV. In the Coulomb breakup, like the radiative capture reaction, the contribution of El transition cross 
section decreases more slowly than E2 cross section with reduction of energy. For instance, at an ad energy of 100 
kcV the the El contribution is less than the E2 one by a factor 70, and at Ead — 10 keV their ratio is equal to 
1/24. This means that the Coulomb breakup cross section is mostly defined by the E2 quadrupole transition at the 
extremely low energy as well. There is a contribution of the E1E2 interference term, but the contribution of the El 
term to the cross section is very small as compared to the E1E2 interference term. 

When we calculate the nuclear breakup cross section according to Eg. pTj) . the radius i?bi of the "black" nuclei ^°®Pb 
is fixed by means of so, that the sum of the Coulomb and nuclear cross sections is close to the experimental data in 
the energy region 0.4 MeV < Ead < 0.73 MeV. For the calculation of we take into account the region where deutron 
emission occurs in the forward direction in the ®Li cm. frame. We do not include the region, Ead < 0.4 MeV, because 
in this region higher-order effects, such as the three-body Coulomb effects [2l|,[23, can give a significant contribution. 
In the region Ead > 0.73 MeV the behavior of the experimental cross section is not clear for the experimental scattering 
angles and the experimental error is comparable with the value of the cross section (see Tables III through VI of Ref. 
[l6|]). The value of ^ minimum when Rti is taken equal to 6.30, 6.94, 6.34, and 6.25 fm corresponding to angles 

©lab = 2°, 3°, 4°, and 6°, respectively. We see that the values of Ru are close to the radius of ^°*Pb. 

The results of calculation of the triple cross sections are shown in Figs. [3] and SI Comparison of our calculated triple 
cross sections with the experimental data shows that the Coulomb breakup gives main contribution into disintegration 
process in the energy range 0.4 MeV < Ead < 0.8 MeV. A discrepancy exists with the experimental data at the energies 
near zero and larger than 0.8 MeV for all scattering angles (particularly in the case 9iab = 3°). As mentioned above 
the purely Coulomb wave function for the ad continuum state maybe applied when relative ad energy is close to 
zero. From analyzing the results done by the authors of Refs. [l^, [Sj, at small values of the Ead relative energy, 
it is clear that the purely Coulomb breaku p c an not explain the behavior of the cross section near the threshold. 
However, the results of the authors of Ref. [l^ showed a huge contribution of the nuclear breakup at all ad energy 
region. Such a phenomenon maybe explained by the wide angular region of the scattering angle of the cm. of ^Li 
(0°-5°) while the grazing angle at this experiment is ^ 2.5°-3°. One expects the purely Rutherford scattering only 
below ~ 2° (see also Figs. 7 and 10 of Ref. (3). At the grazing angle the projectile ®Li moves along the trajectory 
tangential to the surface of the target ^'^^Pb where the nuclear breakup mostly occurs. We note that an analysis of 
the Kiener et al. results concerning angular dependence Ouh at very low Ead energy shows a large divergence with 
the theoretical results following from the purely Coulomb disintegration. The contribution of the nuclear breakup 
increases the total cross section at most by a factor of 2 in the region Ead > 0.8 MeV. A great disagreement still exists 
between the theoretical and experimental results near zero ad relative energy. Figure [S] demonstrates the dependence 
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FIG. 3: (Color online) Panel (a): Energy distribution of the total triple differential cross section for the ^Li breakup as a 
function of the relative energy E^d at an angle of 2° (solid line). The distribution for the purely Coulomb breakup is shown 
with dashed line. The experimental data were taken from Ref. [11]. Panel (6): Same as in (a) but for energy in the close 
vicinity of zero. Negative and positive Ead energies denote backward and forward emission, respectively, of the a particle in 
the "Li cm. frame. Panels (c) and (d): Same as in panels (a) and (6), respectively, but at an scattering angle of 3°. 



of the total cross section and its Coulomb part on the scattering angle at three values of Ead- Figure [5] (a) shows 
that the addition of the nuclear disintegration can not decrease the disagreement with the experimental data at the 
extremely low energies. At still higher energy the inclusion of the nuclear breakup into theory might assist in achieving 
reasonable agreement with the experimental data [Fig. [S](5) and Fig. [5](c)]. Accounting for the higher-order effects 
(as mentioned in I) may lead to better agreement with the experimental data at the extremely low energies. In 
the Kiener et al. work the grazing angle was 13° while the measurement was performed at the scattering angle 
range 2°-6°. Therefore the contribution of the nuclear breakup to the differential cross section was small. The same 
conclusion follows from the results of our calculations presented in Table [J Note that similar behavior is observed for 
other scattering angles. When the energy Ead becomes larger than the resonant energy, the nuclear and Coulomb 
cross sections become comparable to each other for the scattering angles 2° and 3°. However, at scattering angles 4° 
and 6° , the Coulomb breakup cross section remains larger than the nuclear breakup cross section. Figure 3 of Ref. 
[2^ shows how the Coulomb cross section approaches the experimental data when the relative energy Ead increases. 

Figures [3] and m of our paper also depict the existence of the backward- forward asymmetry for deuteron emission in 
the ^Li cm frame especially at smaller scattering angles. Such a type of asymmetry was also observed in the Kiener 
et al. experiment 16]. This asymmetry appears due to the change of sign, both in the Coulomb and nuclear dipole 
transition amplitudes, when the direction of deuteron emission is changed from forward to backward in the ®Li cm 
frame. 

Thus, the results of our calculations confirms that the nuclear contribution exists at the considered scattering angles 
and the ad energies, but it is not as significant as that concluded by the authors of Ref. [l8| . Such a huge contribution 
also contradicts the qualitative estimation given in Ref. (Toj . 
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FIG. 4: (Color online) Same as in Fig. [S] Panels (e) and (/) show results at the scattering angle 4°, while panels (g) and iK) 
show results at the scattering angle 6°. 




FIG. 5: (Color online) Triple differential cross section for the ®Li breakup as a function of the scattering angle 6iab at selected 
values of the relative E^d energy. Solid line is the sum of the Coulomb and nuclear contributions, while dashed line is the 
purely Coulomb one. The experimental data were taken from Ref. pL^ . 



VI. CONCLUSION 



Using the simple two-body approach and the effective potential, we have described the cross section and the 
astrophysical S factor of the a + d — ^Li + 7 radiative capture. The results are in good agreement with the known 
experimental data for the range 0.4 < Ead < 3.0 MeV. For radiative capture the contribution of the El transition 
to the cross section becomes larger than the E2 one at an energy less than 100 keV. The calculated total value of 
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TABLE I: Dependence of the ratio of the nuclear (ctat = fdn"dE ^ ) ^^'^ Coulomb {ac ~ ^cin^dE ^ ) triple cross 
sections for the ^''**Pb(^Li, a d)^''**Pb breakup on Ead energy. The scattering angle Giab is 3°. 



Ead 


(Tjv 






(MeV) 


(mb/MeV/sr^) 


(mb/MeV/sr^) 




1.00 


6.774 


6.316 


1.072 


0.95 


6.372 


5.921 


1.076 


0.90 


5.932 


5.376 


1.104 


0.85 


5.459 


4.641 


1.176 


0.80 


4.981 


3.809 


1.308 


0.77 


4.806 


3.968 


1.211 


0.76 


4.847 


4.708 


1.030 


0.75 


5.048 


6.632 


0.761 


0.74 


5.663 


11.73 


0.483 


0.73 


7.555 


26.98 


0.280 


0.72 


13.97 


79.58 


0.175 


0.71 


24.61 


172.8 


0.142 


0.70 


16.56 


111.9 


0.148 


0.65 


5.043 


18.10 


0.279 


0.60 


3.757 


10.49 


0.358 


0.55 


3.007 


7.621 


0.394 


0.50 


2.395 


5.854 


0.409 


0.45 


1.856 


4.503 


0.412 


0.40 


1.378 


3.372 


0.409 


0.35 


0.964 


2.400 


0.402 


0.30 


0.621 


1.581 


0.393 


0.25 


0.354 


0.927 


0.382 


0.20 


0.168 


0.453 


0.371 


0.15 


0.058 


0.162 


0.360 


0.10 


1.07x10"^ 


3.060x10"^ 


0.349 


0.05 


3.05x10"* 


9.047x10"'' 


0.337 



the astrophysical S factor equals to S{0) = 1.22 MeVnb, while Sei{0) = 1.01 MeVnb (83% of the total S factor) 
and Se2{0) = 0.21 MeVnb (17% of the total S factor). The results for the purely Coulomb breakup are in good 
agreement with the known experimental data for the range 0.4 < Ead < 0.8 MeV. This shows the validity of the 
semiclassical method of the calculation for the Coulomb breakup at high energy collision. The contribution of the El 
transition in the Coulomb breakup is always less than the E2 one for all energy regions. The nuclear disintegration is 
analyzed by diffraction method which can be applied at small scattering angles. The radius of "black" target ^osp]-, 
is taken as a fit parameter and the application of the method gives reasonable fit to the radius of the target. 
Our calculated nuclear distortion is not large and suggests an overestimation of the contribution of nuclear distortion 



made by Hammache et al. jl8| . A comparison of our calculation of the triple cross section, consisting of the Coulomb 



and nuclear parts, with the experimental data, shows that disagreement still exists for Ead near zero. Taking into 
account the higher-order effects may reduce this discrepancy at the low ad relative energy {Ead < 0.3 MeV). For 
instance, the three-body Coulomb effects are known to strengthen with decreasing ad relative energy [21I, [S^ . Due to 
the existence of the nuclear breakup and the small contribution of the El Coulomb disintegration to the total cross 
section, it is impossible to extract the correct value of the astrophysical S factor for a + d ^ ^Li-l-7 radiative capture 
at low Ead energy. Extrapolation to zero energy can give the value of the E2 component, if one is able to separate the 
contributions of the nuclear disintegration and the higher-order effects from experimental data. To us, the simplest 
way to get an accurate value of the astrophysical S factor of peripheral reactions, is to measure the ANC of the wave 
function in the bound state with a high accuracy which governs the overall normalization of the peripheral reaction 
cross section near zero energy. 
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